Quantization of perturbations during inflation in the 1 + 3 covariant formalism 
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This note derives the analogue of the Mukhanov-Sasaki variables both for scalar and tensor 
perturbations in the 1+3 covariant formalism. The possibility of generalizing them to non-flat 
Friedmann-Lemaitre universes is discussed. 
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The theory of cosmological perturbations is a corner- 
stone of the modern cosmological model. It contains two 
distinct features. First, it describes the growth of den- 
sity perturbations from an initial state in an expanding 
universe filled with matter and radiation, taking into ac- 
count the evolution from the super-Hubble to the sub- 
Hubble regime. Second, it finds the origin of the initial 
power spectrum of density perturbations in the amplifi- 
cation of vacuum quantum fluctuations of a scalar field 
during inflation. 

Two formalisms are used to describe the evolution of 
perturbations. The first relies on the parametrization 
of the most general spacetime close to a Friedmann- 
Lemaitre (FL) universe and on the construction of gauge 
invariant variables The second is based on ageneral 
1 + 3 decomposition of the Einstein equation |2[]. The 
philosophies and advantages of these two formalisms are 
different. 

The Bardeen formalism is restricted to perturbations 
around a FL universe (with metric g^ v ) and expands the 
spacetime metric as g^ v = g^ v +7^. The metric pertur- 
bation is then decomposed as 

lp , v A^Ax v = 2a 2 (77) [-Adrj 2 + [D t B + B % ) dsMr? 

+ (C 7 y + DiDjE + D (l Ej) + Ekj) ditW] (1) 

and the 10 degrees of freedom are decomposed in 4 scalars 
(A, C, B, E), 4 vectors {E,,B t with D t B l = 0, . . .) and 2 
tensors (.Ey with E\ = D^E^ — 0). During inflation, it 
was shown that the Mukhanov-Sasaki (MS) variable v, 
a gauge invariant variable that mixes matter and metric 
perturbations, must be quantized [1,0, Q. This approach 
is thus completely predictive (initial conditions and per- 
turbations evolution) for an almost flat FL spacetime. 
It is not straightforward to extend it to less symmetrical 
inflationary models @, 0] , for which the quantizationpro- 
cedure has not been investigated (see however Ref. [6[ for 
non flat FL universes and Ref. [?J for non-FL spacetimes) . 

The 1 + 3 covariant description assumes the existence 
of a preferred congruence of worldlines representing the 
average motion of matter. The central object is the 4- 
velocity u a of these worldlines together with the kine- 



matical quantities arising from the decomposition 



V a Mfc = -u a u b 



-Qh a b + (T a b + UJab (2) 



where h a b = g a b + u a u b- © = V a u a is the expansion rate, 
a a b the shear (symmetric trace-free with a a bU a = 0), 
uj a b the vorticity (antisymmetric with uj a bU a — 0) and 
u a = u b \7 b u a ). A fully orthogonally projected covariant 
derivative D a , referred to as a spatial gradient, is defined 
and a complete set of evolution equations can be obtained 
for these quantities (together with the electric and mag- 
netic parts of the Weyl tensor, E a b and H a b and the mat- 
ter variables) without needing to specify the spacetime 
geometry (see e.g. Ref. [lfj for details). The formalism 
can then be used to study the evolution of perturbations 
in various spacetimes (see e.g. Ref. [TT| for Bianchi uni- 
verses) and be yon d the linear order in the case of almost 
FL universes [12|, O)- ^ n this case, the perturbation 
variables have a clear interpretation and have been re- 
lated to the Bardeen variables both for fluids [3] and 
scalar fields On the other hand, the scalar-vector- 

tensor decomposition is not straightforward [l6| and the 
analogue of the Mukhanov-Sasaki variable v has not been 
derived so that it is difficult to argue which quantity must 
be quantized in this formalism. 

The goal of this article is to identify this variable in the 
1 + 3 covariant formalism. In § |T1 we recall the construc- 
tion of the Mukhanov-Sasaki variable. We then identify 
in § [TT] its counterpart in the 1 + 3 covariant formalism 
for an almost FL spacetime. § IIIII adresses the question 
of gravitational waves. Finally we conclude on the use 
and possible extensions of these variables in § IIVI 



I. QUANTIZATION OF THE 
MUKHANOV-SASAKI VARIABLE 

Focusing on scalar modes of the decomposition |T]), 
one defines the two gauge invariant potentials, <& = A + 
H(B - E') + (B — E'Y and * = -C - H(B - E') where 
H = a! I a. For a universe filled by a single scalar field 
if, one can define the gauge invariant field fluctuation 
Q = Sip — (p'C/H. The curvature perturbation in the 
comoving gauge is then given by 
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when assuming a background FL universe with flat spa- 
tial sections. Introducing the two variables 



2a$ 

Kip' 



H 

cup' 



H 

aip 1 



(4) 



with k = 8irG, 1Z takes the simple form 1Z = 9u' — Q'u, 
where use has been made of the background equation 
9' = (-aQ<p" +W- H 2 ) /atp' = -Op"/p' - up '/2a. If 
we define 



zTZ = aQ 



(5) 



it reduces to the simple relation zv = (uz)'. Thus, TV 
takes the simple form TV = 9u" — 9"u. Additionnaly, 
the perturbed Einstein equations imply that TV = 9Au. 
Provided ip' + 1 0, i.e. we are not in a strictly de Sitter 
phase, this leads to 







which can also be written as 

z 2 TZ' — All TZz 2 dr] 



(6) 



(7) 



This form is strictly equivalent to K" - ATl+2(z' / z)W = 
0, that is to 



v" - ( A + — ] v = 0. 



(8) 



This is the equation of a harmonic oscillator with time 
varying mass m 2 — z" / z. To show that v is indeed the 
canonical variable to be quantized (contrary to u which 
satisfies a similar equation with z replaced by 9), one 
needs to expand the action of the scalar field coupled to 
gravity up to second order ||. The standard procedure 
is then to promote v to the status of an operator and set 
the initial conditions by requiring that it is in its vac- 
uum state on sub-Hubble scales (in Fourier space, that is 
when krj — * -co). At the end of inflation, all cosmologi- 
cally relevant scales are super-Hubble {krj <C 1) and the 
conservation of 1Z on these scales is used to propagate 
the large scale perturbations generated during inflation 
to the post-inflationary eras. 



II. QUANTIZATION IN THE 1+3 FORMALISM 

A. Generalities 

Let us first stress an important point concerning time 
derivative. In the coordinate based approach, the dot 
refers to a derivative with respect to the cosmic time. 
In the covariant formalism, a natural time derivative 
is introduced as u a V a which is a derivative along the 
worldline of the observer. For an almost FL space- 
time, it is clear that for any order 1 scalar quantity X, 



u a V a X — dt X since u a has to be evaluated at the back- 
ground level. It is indeed not the case anymore for vec- 
tor or tensor quantities, or at second order in the per- 
turbations. Another time derivative can be constructed 
from the Lie derivative along u a , C u . This derivative 
exactly matches the derivative with respect to cosmic 
time for any type of perturbations [15]. For instance, 
C u X a — u b VbX a +XbV a u b , and it is easily checked that 
the second term exactly compensates the term arising 
from the divergence of u a . We recall that the Lie deriva- 
tive satisfies 



D a {X) = C u (D a X) — u a X. 



(9) 



In the 1 + 3 covariant formalism, the main equations 
for a spacetime with vanishing vorticity {uj a b — in the 
decomposition ([2])) are the Raychaudhuri equation 







1 



9 2 = D a u a - -( p + 3P) -2a 2 



3 "2 
the Gauss-Codacci equation 
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together with the conservation equations 

p + 9(p + P) = 0, D a P=-{p + P)u a 



(10) 



(11) 



(12) 



where 2a 2 = a a b<J ab . For a scalar field p = ip 2 /2 + V and 
P = t/j 2 /2 — V with ip = <j>. We define the scale factor S 
by the relation S / S = H = 0/3. The previous equations 
imply that 



H 



i 



2 T 6 
and the Gauss equation leads to 



K + D a ii a 



t = --G(/C + 2D a u a ) + 2(cr 2 )- + 26cr 2 



and 



D a IC 



-6D a + 2KipD a il> + 2D a (a 2 ). 



(13) 



(14) 



(15) 



Let us now introduce the central quantity in our dis- 
cussion 



7^4 



J (D a O)dT-D a (J 6d 



^ (16) 



where r is the proper time along the fluid flow lines. It 
follows from the identity © that 



£"uTZ a 



-Hii n 



H(D a ^)/^, 



(17) 



which is a useful relation in order to close the equa- 
tions in TZ a . The conservation equation implies that 



3 



£ u (ipD a ilj) = —ipD a (<dip) which, combined with Eq. (fT5f . 
gives 



l S 2 D a (IC-2a 2 ) = Ua + ^ 

4 6K1p Z 



CuU a + ^QU a 



(18) 



with U a = nS 2, 4iD a 'ip /2. Now, introducing 
V a = zS 2 



1 D a (IC~2a 2 ) + (* 2 - 1 -IC-D b u b )^. 

4 b ip 



(19) 



and developping C u {zU a ) with the help of Eqs. |T 
and (1131). it can be shown that 



1 

Sz- 



zV a 



dr. 



(20) 



Whatever the values taken by the shear and the spatial 
curvature, this intermediate result is exact, that is valid 
for any order of perturbation. Note its similarity with 
Eq. ©. 



S 2 D a KL 



S 2 D b D b K a = ATZ a 



(22) 



up to a constant F a satisfying C u F a = which can be 
absorbed in the integration initial boundary surface of 
Eq. ([16]). Thus Eq. d20j) reads 



l 

~sl 2 



-ATZ a dr. 



(23) 



At this stage, it is useful to introduce a vector field w a = 
Su a , and the conformal proper time f defined by Sdr — 
dr. It is easily seen that for a spatial vector (i.e. u a X a — 
0) C w X a = SC u X a . The Lie derivative along w a matches 
at first order the derivative with respect to the conformal 
time 77, as the Lie derivative along u a was matching the 
derivative with respect to the cosmic time. For scalars we 
thus use the notation X' = C W X. With this definition, 
Eq. (j2"3"|) can be recast as 



B. Flat FL spacetimes 

Let us now focus on homogeneous spacetimes with flat 
spatial sections; this includes the FL spacetimes for which 
the shear vanishes and Bianchi I spacetimes. Homegene- 
ity implies that the spatial gradient of any scalar func- 
tion vanishes (D a f = 0, V/) and in particular D a P = 
so that ii a = 0. Flatness implies that JC — and that the 
3D-Ricci tensor vanishes, ( - 3 - > R a b — 0, which leads to the 
simplified equation for the shear 

(Tab + &°ab = => {a 2 ) ■ + Q<7 2 = 0. (21) 

This means that the only non-vanishing quantities are 
0, (f>, ip and cr. At first order in the perturbations, we 
have to consider the gradients of these quantities but also 
terms like u a and K,. 

FL spacetimes are also isotropic. This implies that a = 
and S = a for the background. Consequently we can 
discard gradients of a 2 which are second order. Thus, 
the only remaining term in the expression (|19[) of V a is 
\zS 2 D a K. 

In order to get a closed equation for lZ a , we need to 
express V a in terms of lZ a in Eq. (|20| . Taking the spatial 
gradient of Eq. (fT4| and using Eq. (J9j) for handling the 
time derivative, we obtain the first order relation 

where the last equality follows from the flat background 
assumption. Eq. (|17p and the commutation relation 
S 2 D b D b C u X a = C u (S 2 D b D b X a ), valid at first order, 

imply C u ( s2] l alC ^ = C u (S 2 D b D b K a ) which once inte- 
grated leads to 



z 2 C w TZ a = J z 2 AK a dT, (24) 

which is similar to Eq. ([7]). By the same token, we de- 
duce that v a defined in Eq. (TTS)) satisfies Eq. ([5]). It can 
be checked that its spatial components are linked to the 
MS variable at first order in perturbations by vi = diV, 
and consequently the initial conditions obtained from the 
quantization of v can be used to set the initial conditions 
for v a and then lZ a . v a is thus the analogue in the 1 + 3 
formalism of the MS variable v in the Bardeen formalism 
and it satisfies 

Cl (O- (A+— )v a = 0, (25) 



III. GRAVITATIONAL WAVES 

It can be shown that the magnetic part of the Weyl 
tensor H ab is a good variable to describe the gravitational 
waves [ill [3], and it satisfies at first order for flat FL 
spacetimes 

C 2 w H ab + 2C W (HH ab ) ~ AH ab = . (26) 

£ a b defined by C w £ a b = H a b satisfies 

C? w £ab + 2H£ W (Sab) - A£ab = (27) 

where the integration constant is set to as for the scalar 
case. This is exactly the equation satisfied by the grav- 
itational waves Eij. Indeed, these variables are linked 
at first order by £ a b — e cd < a d c E b> d [3, where e ab c is 
a completely antisymmetric tensor normalized such that 



4 



e 12 3 = 1. With Ma6 
order to 



;=£ a b, Eq. (HH) leads at first 



S" 



C w {/J,ab) - ( A + — ) /i ah = 



(28) 



which is the equation for an harmonic oscillator with a 
time varying mass S" / S. Just as for scalar perturbations 
case, the quantization of the gravitational waves in per- 
turbation theory, can be used to set the initial conditions 
for ji a b and thus for £ a b and H a b- 



IV. CONCLUSIONS 

We have identified the scalar and tensor variables 
that map to the Mukhanov-Sasaki variables when con- 
sidering an almost-FL universe with Euclidean spatial 
sections in the 1 + 3 covariant formalism. Let us 
stress that in Ref. [Hj], H a = ~D a a (where a = 
\ J c Qdr] is the integrated volume expansion along u a ) 
was proposed. But clearly, this maps at linear order to 
d a (Tl - 3 / A(V + E')dij) , where V is the scalar part of 
the perturbation of the velocity. The additional term 
in our definition (|16[) cancels this discrepancy as it can 
be seen from the constraint |-D a 9 = D a b a which im- 
plies that at first order D a Q = ^d a AV. Alternatively, 
this can be seen directly on the expression of in terms 
of perturbation variables by use of the (0 — i) Einstein 
equation. The two variables agree at leading order on 
super-Hubble scales. 

Two generalizations with less restrictive backgrounds 
can be considered: flat but anisotropic spatial sections 



(K = = 0, cr ^ 0), and isotropic but non-flat FL 
spacetimes (K ^ 0, a = 0). The cornerstone of the 
derivation of § [HI is the possibility of expressing V a only 



in terms of 7Z a to get a closed equation from Eq. (|20|) . 
In the first case, two other terms in Eq. [fT9|) contribute 
at first order. zS 2 a 2 ^^- changes the definition of the 

effective varying mass, but the term ^-D a {a 2 ) acts as 
a source in the R.H.S of Eq. ([5]). This term represents 
a coupling of the gravitational waves (a a b at first order) 
with the shear of the background spacetime. As for £ ab , 
Eq. (|27|) will be supplemented with an integral non trivial 
source term which couples the b ackg round shear to the 
electric part of the Weyl tensor pj|. Thus both equa- 
tions are mixed with these new source terms, which are 
of the same order of magnitude as the quantized vari- 
ables. Note that this not surprising since at second order 
in the perturbations around a FL spacetime the scalar 
and tensor degrees of freedom are coupled [H, [2(| • 
In the second case (K ^ 0, a = 0), 



V a 



-D a IC - K 



-C 



(29) 



The spatial gradient of Eq. (|14p leads at first order to 



C-uC a — 
where £ a = D a a 



S 2 H 



D a (D b D b i;) - 3K£ u ( a 



(30) 



is a possible nonlinear general- 



ization of the curvature perturbation on uniform density 
hypersurfaces. From the conservation equation (|12[) it 
can be shown [l3|, [TtJ that £ a is conserved (in the sense 
of the Lie derivative) on large scales for adiabatic pertur- 
bations, and thus C a is also conserved in the same sense 
on super-Hubble scales for adiabatic perturbations 
Because of this term, C a cannot be expressed solely in 
terms of lZ a as it has been done with Eq. ()22l) . Indeed, 
it contains an additional term involving £ a [12J. 
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